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1.  Objectives 


The  objectives  of  the  research  were  as  follows. 

(1)  Develop  an  accurate  and  general  framework  of  smart  composite  structures  including 

•  Accurate  representation  of  through  the  thickness  shear  deformation. 

•  Presence  of  multiple  delaminations. 

•  Two-way  coupling  between  electrical,  mechanical  and  thermal  fields  in  the  constitutive 
relations 

(2)  Investigate  the  use  of  dynamic  strain  and  RMS  values  in  damage  detection.  Design  a  control 
system  to  minimize  the  effect  of  damage. 

2.  Accomplishments 

A  general  framework  has  been  developed  for  the  analysis  of  laminated  smart  structures 
with  embedded  and/or  surface  bonded  piezoelectric  actuators  and  sensors.  The  theory  addresses 
the  following  issues:  (i)  accurate  description  of  transverse  shear  effects  in  the  laminated 
structural  model,  (ii)  presence  of  multiple  delaminations  and  (iii)  characterization  of 
delamination  placement  and  size  on  dynamic  characteristics.  The  mathematical  model  has  been 
implemented  using  finite  element  method  and  extensive  numerical  investigations  have  been 
performed  to  validate  the  model  with  three-dimensional  results  and  available  experimental  data. 
Results  indicate  that  the  presence  of  delaminations  affect  the  modal  strain  distributions  in  smart 
composite  laminates.  The  developed  theory  has  been  used  to  characterize  delaminations  in  smart 
composite  plates  of  arbitrary  thickness  and  results  indicate  that  dynamic  strain  is  a  better 
indicator  of  damage  than  frequencies  or  mode  shapes.  Significant  differences  are  noted  in  the 
RMS  values  of  the  response,  calculated  using  modal  strain,  indicating  the  presence  of  damage. 
A  control  system  has  also  been  designed  to  minimize  the  effect  of  damage.  New  damage  indices 
based  on  modal  strain  are  proposed  and  show  promising  results. 
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3.  Development  of  a  Higher-Order  Based  Finite  Element  Model 


Fig.  1  Laminate  cross  section  with  delamination 

The  general  higher  order  displacement  field  has  been  extended  to  model  composite  plates 
with  surface  bonded  piezoelectric  materials  and  delaminations  (Fig.  1).  The  inplane 
displacements  are  effectively  expressed  by  a  cubic  function  through  the  thickness  (z)  and  the 
transverse  displacement  is  assumed  to  be  independent  of  z.  To  model  delamination  in  such 
structures,  it  is  necessary  to  partition  the  laminate  into  several  different  regions  as  shown  in  Fig. 
1.  These  regions  include  the  undelaminated  region,  the  region  above  the  delamination,  and  the 
region  below  the  delamination.  The  interface  between  the  undelaminated  region  and  the 
delaminated  region,  indicated  by  the  dashed  line  in  Fig.  1,  is  denoted  S.  The  general  form  of  the 
higher  order  displacement  field  is  independently  applied  to  each  of  these  regions  to  describe 
displacements,  which  account  for  slipping  and  separation  due  to  the  delamination.  The  higher 
order  displacement  theory  is  used  to  model  a  composite  laminate  with  multiple  delaminations. 

u,  (x,  y,  z, l)  =  u0i  (x,  y, l)  +  (z-  zcl  )[-w0l>  +  a f  (x, y,  /)]  +  (z  -  zd  )2  uv (x,  y,  t)  +  (r  -  zci )3  u3i  (x,  y,  t) 

v,  (x,  y,  z, t)  =  v0i  (x,  y, /)  +  (z  -  zd  )[-»%,,  +  /?,  (x,  y,  <)]  +  (z  -  zci  )2  v2(  (x,  y, t)  +  (z-  zci  )3  vJ;  (x,  y,  t)  (1) 

wi(x,y,z,t)  =  wai(x,y,t) 

where  m,-,  v,-  and  w,  are  the  displacement  functions,  uo\,  vo,  and  wm  denote  the  midplane 
displacements  of  a  point  (x,y),  wm,x  and  woj,y  represent  the  rotation  of  normals  to  the  midplanes 

and  a,  and  are  the  additional  rotations  due  to  shear  deformation.  The  quantities  uu,  V2i  and 

V3i  represent  higher  order  functions.  The  index  i  is  used  to  designate  the  three  regions,  i-1  the 
undelaminated  region,  i-2  the  region  above  delamination  and  i=3  the  region  below 
delamination,  while  zci  are  the  positions  of  the  midplane  for  each  region  with  respect  to  laminate 
midplane  (Fig.  1)  and  t  denotes  time. 
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The  transverse  shear  stresses,  ov  and  05,  must  vanish  at  all  free  surfaces,  including  at 
delamination  interfaces.  This  leads  to  the  following  refined  displacement  field. 

u,  (x,  y,  z,  0  =  u0l  (*,  y,t)  +  (z  -  zd  )[-  v% ,  +  a,  (x,  y,t)]-4(z- zd  ) 3  a,  (x,  y,  t)K}hf ) 

v, .  (x,  y, z,  t )  =  v0i  O,  y,  t)  +  (z  -  zd  )[-w0i  y  +  (i,  (x,  y,t) ]-  4(r  -  zc()3  /?,  ( x,y,t)/(3h f  )  (2) 

wl(x,y,z,t)  =  w0l(x,y,t) 

In  addition,  continuity  conditions  of  displacement  and  traction  are  imposed  at  the 
midplane  points  of  each  region  on  delamination  lateral  boundaries  S  (Fig.  1)  as  follows. 

u, (x,y,  zd ,  0  =  u,  (x,  y,  zd ,  t) 

v,  (x,y,zd,t)  =  vi(x,y,zci,t)  '  =  2,3  (3) 

wt  (x,y,zcl,t)  =  w,(x,y,zci,t) 

and 

0„i(x,y>zcl,t)  =  <rni(x,y,zd,t) 

r„/,i  (x,y,  zc2 ,  t)  =  tm(x,  y,  zd ,  t)  i  =  2,3  (4) 

r„vi  (*>  y,  zc2 , 0  =  rnvl  (x,y,  zd ,  /) 

In  the  finite  element  formulation  these  continuity  conditions  are  to  be  satisfied  at  the 
mesh  node  locations.  The  constitutive  equations  take  the  following  vector-matrix  form. 

i  =  2,3  (5) 

where  Us  represents  the  generalized  displacement  vector  in  the  regions  above  and  below 
delamination,  Us  is  the  displacement  vector  in  the  undelamianted  region  and  T  is  the 
transformation  matrix  between  regions. 

In  the  context  of  laminated  composite  plates,  the  constitutive  equations  for  the  k!h 
laminae  take  the  following  form. 

c*  =  Qkek  (6) 

where 

o-k=k  ^  (?) 

*k=[Ex  £y  Y  yz  Yxz  YxyJ  ^ 

and 
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(9) 


en  e>2  o  o  e,*' 

2l2  022  ®  ®  C?26 

2a  =  0  0  a4  245  0 

0  o  245  2s5  0 

_2l6  226  ®  ®  2 66 _  I 

In  the  above  equations,  <x  and  e  denote  the  stress  and  the  mechanical  strain  in  the  material 
coordinates  respectively,  Q  is  the  elastic  stiffness  matrix. 

The  finite  element  method  is  used  to  implement  the  theory.  The  equations  of  motion  are 
derived  using  Hamilton’s  principle  in  the  discrete  form  and  can  be  written  as  follows. 

Mij  +  Kq  =  F  (10) 

where  M,  K,  F  and  q  denote  the  mass  matrix,  the  stiffness  matrix,  the  force  vector  due  to  a 
distributed  load  and  the  nodal  displacement  vector,  respectively.  Bilinear  shape  functions  are 
used  for  interpolating  the  inplane  displacements  u  and  v  and  rotations  a  and  while  a  16  term 
cubic  polynomial  is  used  for  the  transverse  displacement  w.  The  representative  four-noded 
element  has  32  degrees  of  freedom. 

Damage  Detection  and  Vibration  Control 

To  characterize  the  dynamic  characteristics  of  smart  composite  plates,  both  delaminated 
and  nondelaminated,  are  investigated  next.  The  plate  is  subject  to  impulse  and  RMS  values  of 
the  response  are  calculated  in  an  attempt  to  understand  the  mechanics  of  damage  detection.  The 
equations  of  motion  can  be  expressed  as  the  following  state  equations. 

{^  =  [A]{x}  +  [B]{u}  +  {Bw}w 
{y}  =  [C]{x} 

where  [A],  [B]  and  [C]  are  the  state  matrix,  the  control  matrix  and  the  output  matrix  respectively. 
The  quantity  w  is  white  noise  and  {Bw}  is  noise  vector.  After  vibration  analysis,  a  modal 
reduction  is  performed  using  the  first  ten  elastic  modes.  Therefore,  the  state  matrix  is  of  size 
(20x20).  The  RMS  values  of  response  for  the  elastic  modes  are  obtained  from  the  following 
equation. 

or?  =[[C][X][C]T]H,i  =  l,2,..,Ns  (12) 
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where  [C]  is  the  output  matrix,  [X]  is  the  state  covariance  matrix  and  Ns  represents  the  number 
of  structural  states.  The  state  covariance  matrix  [X]  of  the  system  is  the  solution  of  a  Lyapunov 
equation  of  the  form 

[A][X]  +  [X][Af  +  {Bw}Qw{Bw}t  =0  (13) 

where  Qw  is  the  intensity  of  white  noise.  The  RMS  values  of  the  elastic  modes  can  be  converted  into 
displacement  and  velocity,  at  the  nodal  points,  by  multiplying  the  covariance  matrix  by  the  matrix 
<t>2  as  follows 

T  TO  0] 

=  [%*■*.  =[0  <14) 

where  aq(y  represents  the  RMS  values  of  nodal  displacement  and  nodal  velocity  and  <t>  is  the 

eigenmatrix.  The  RMS  values  of  the  nodal  displacements  due  to  the  disturbance  are  calculated 
for  both  the  undelaminated  and  the  delaminated  plate  models. 

Active  Control  System  Design  to  Minimize  the  Effect  of  Damage 

The  damage  detection  algorithm  is  extended  to  design  an  active  control  system.  The  goal 
is  to  reduce  the  RMS  values  and  thereby  minimize  the  effect  of  such  damage.  It  is  assumed  that 
all  states  are  available  for  feedback.  The  gain  matrix  of  the  completely  controllable  and 
observable  linear  dynamic  system  is  expressed  as  follows 

{u}  =  -[KG]{x}  (15) 

The  right  and  left  eigenvalue  problems  of  the  closed  loop  system  can  be  written  as 

([A]-[B][KG]){f },  =  Wo  ([A]-[B][KG])T{M/C}i  =  (16) 

where  {<J>c}i  and  are  the  right  and  left  eigenvectors  of  the  closed  loop  system,  respectively, 

corresponding  to  the  eigenvalue  . 

The  pole  placement  algorithm,  based  on  Sylvester  equation,  utilizes  the  parameter  vector 
defined  by 


{h}  j  =  [KG]{<t>c}j 


(17) 
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Substitution  of  Eq.  (17)  into  Eq.  (16)  yields  the  following  Sylvester  equation 


([Al-QlDttVtBlfi},  (IB) 

or  in  matrix  form 

[A][OJ-[Oc][Ad]  =  [B][H]  (19) 

where 

[Ad]  =  Diag(A,c1,  A,c2,  (20) 

[H]  =  [{h}[,  {h}2, ...,{h}n] 


and  [Ad]  contains  the  desired  closed  loop  eigenvalues.  From  Eq.  (18),  it  can  be  seen  that 

{f}i  =  ([A]-^[I])'1[B]{h}i  (21) 

If  X]  are  distinct  from  their  open  loop  positions,  the  columns  of  [H]  generate  the  corresponding 
closed  loop  eigenvectors.  For  the  given  parameter  matrix  [H],  the  closed  loop  modal  matrix 
[<J>  J  can  be  determined  by  solving  Sylvester  equation.  The  gain  matrix  can  be  obtained  from 

Eq.  (17)  as  follows 

[KG]=[H][a>crl  (22) 


Development  of  Damage  Index 

Finally,  a  new  damage  index  is  formulated  using  the  difference  between  the  inplane 
modal  strains  of  the  delaminated  structure  and  those  of  a  similar  healthy  structure.  The 
difference  in  modal  strains  is  squared  to  diminish  the  effects  of  numerical  errors  or  experimental 
noise.  This  procedure  requires  an  undamaged  reference  but  is  valid  for  any  type  of  structure 
including  smart  laminate  with  stiffness  discontinuity.  The  new  damage  index  is  defined  as 
follows. 

5m,=k-(4)J2  (23> 
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where  ml  is  the  number  of  the  element  under  investigation,  r  is  the  number  of  vibration  mode  for 
which  the  modal  strain  is  computed  and  k=  1,2,6  is  the  appropriate  modal  strain. 

Results  And  Discussions 

Results  are  presented  for  Graphite/Epoxy  composite  plates  with  five  pairs  of  surface 
bonded  piezoelectric  actuators.  The  plates  are  clamped  along  one  edge  and  other  edges  are  free. 
The  dynamic  responses  are  investigated  for  plates  with  and  without  a  small  delamination  (Fig. 
3).  The  stacking  sequence  is  [0°/0°/90°]s  and  each  ply  has  a  uniform  thickness  of  0.0000134m. 
The  plate  length  is  0.305m,  width  is  0.076m  and  total  thickness  is  0.000804m.  A  delamination  is 
placed  at  midplane  in  the  region  between  0.0152<x<0.0456  m,  0.0305<y<0.061  m.  The 
material  properties  are:  Ei=  98.0  GPa,  Ej-  7.9  GPa,  vn=  0.28,  G/2=  Gu=  5.6  GPa,  G23-  2.4 
GPa,  p=  1520  Kg/m3  for  the  composite  plate  and  E=  63  GPa,  v=  0.31,  G=  24.2  GPa,  p=  5000 
Kg/m3,  di2=  250xl0'n  m/V  for  the  piezoelectric  material.  Table  1  shows  the  effect  of 
delamination  on  the  first  six  natural  frequencies  of  the  composite  plate.  Some  changes  are 
observed  in  the  values  of  the  first  three  frequencies  as  a  result  of  the  small  reduction  in  the 
structural  stiffness  due  to  delamination. 


Fig.  2  Delaminated  composite  plate  with  piezoelectric  actuators. 

The  RMS  values  of  the  nodal  displacements  due  to  the  disturbance  are  calculated  for  both 
nondelaminated  and  delaminated  plates.  Figures  3  and  4  show  the  responses  of  the  inplane 
displacement  u,  due  to  disturbance,  for  the  nondelaminated  and  the  delaminated  plates,  respectively. 
As  seen  from  these  figures,  very  significant  differences  are  observed  in  the  RMS  values  of  inplane 
displacement  between  cases  without  and  with  delamination.  In  the  presence  of  delamination,  the 
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RMS  values  of  the  response  due  to  disturbance  increases,  as  shown  in  Fig.  4  compared  to  the 
nondelaminated  plate  (Fig.  3).  This  is  once  again  due  to  the  reduction  in  plate  stiffness  caused  by 
the  presence  of  delamination.  Significant  jumps  in  the  RMS  values  in  the  delamination  boundaries 
are  also  observed.  This  information  can  be  used  to  identify  delmaination  location  and  size.  The 
RMS  values  of  the  out  of  plane  displacement  (w)  were  also  calculated.  However,  it  was  found  that 
the  changes  in  RMS  distributions  of  this  displacement,  between  the  delaminated  and  the 
nondelaminated  plates,  were  not  significant  enough.  Therefore,  this  information  cannot  be  used  to 
accurately  capture  the  delamination  zone. 


Nondel. 

Plate  (Hz) 

Del. 

Plate  (Hz) 

1 

11.7 

11.2 

2 

63.0 

62.4 

3 

90.1 

86.7 

4 

174 

173 

5 

196 

190 

6 

331 

326 

Table  1  Natural  frequencies  of  nondelaminated  and  delaminated  plates. 


Fig.  3  RMS  values  of  inplane  displacement  ,u,  of  nondelaminated  plate. 
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Fig.  4  RMS  values  of  inplane  displacement,  u,  of  the  delaminated  plate. 


Nondelaminated  plate 

Delaminated  plate 

Delaminated  plate 

(open  loop) 

(open  loop) 

(closed  loop) 

-.1475+73.771 

-.1410170.511 

-7.377173.771 

-.7921+396.71 

-.78431392.11 

-39.601396.01 

-1.1311565.91 

-1.089+544.81 

-56.59+565.91 

-2.19011095.1 

-2.17611088.1 

-109.511095.1 

-2.46311231.1 

-2.39011195.1 

-123.111231.1 

-4.16212081.1 

-4.10612053.1 

-208.112081.1 

-4.986+2493.1 

-4.90612453.1 

-249.312493.1 

-5.67312836.1 

-5.60612803.1 

-283.612836.1 

-6.357+3178.1 

-6.094+3047.1 

-317.813178.1 

-6.39513197.1 

-6.38013190.1 

-319.713197.1 

Table  2  Changes  in  eigenvalues  of  the  open  loop  and  closed  loop  systems. 
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Table  2  presents  the  changes  in  eigenvalues  for  both  delaminated  and  nondelaminated  plates 
with  the  application  of  the  active  control  system.  As  seen  from  this  Table,  significant  increases  are 
observed  in  the  real  part  of  the  closed  loop  eigenvalues  for  the  delaminated  plate.  This  indicates 
improved  damping  with  closed  loop  control.  The  imaginary  parts  of  the  closed  loop  eigenvalues  for 
the  delaminated  plate  are  almost  identical  to  the  corresponding  open  loop  values  for  the 
nondelaminated  plate.  This  implies  that  the  natural  frequencies  of  the  two  plates  are  almost 
identical.  The  RMS  values  of  the  first  three  elastic  modes  corresponding  to  the  open  loop  system  of 
the  nondelaminated  and  delaminated  plates  are  (15.227,  0.0292,  0.2093)  and  (21.035,  0.0142, 
0.2578),  respectively.  With  the  application  of  the  active  control  system,  these  values,  for  the 
delaminated  plate,  change  to  (0.681,  0.091,  0.0051).  This  indicates  a  very  significant  reduction  (up 
to  78  percent)  in  the  magnitude  of  the  RMS  values  for  the  first  and  third  modes.  The  RMS  value  of 
the  second  mode  increases  compared  to  the  open  loop  case.  However,  the  contribution  of  this  mode 
is  not  significant  enough  to  cause  stability  problem  or  degradation  in  the  performance  of  the  active 
control  system.  Therefore,  the  controller  is  successful  in  improving  the  dynamic  characteristics  of 
the  delaminated  plate. 

To  investigate  the  control  performance  of  the  designed  system,  the  time  histories  of  the  tip 
displacement  and  the  control  voltages  applied  to  the  actuators  are  calculated.  Figure  5  shows  the 
time  histories  of  the  open  loop  and  closed  loop  systems,  due  to  disturbances,  for  the  nondelaminated 
and  the  delaminated  plates.  The  control  system  significantly  improves  the  response  of  the 
delaminated  plate  compared  to  the  nondelaminated  plate  (open  loop).  Figure  6  presents  the  control 
voltage  applied  to  the  actuators  due  to  the  disturbance.  It  can  be  seen  that  maximum  applied  voltage 
is  less  than  200  Volts  implying  no  saturation  effects. 

Next,  composite/smart  composite  plates  with  single  or  multiple  delaminations  are 
analyzed  to  verify  the  performances  of  the  proposed  damage  indices.  The  geometry  of  the  test 
structure  consists  of  eight  plies,  each  with  thickness  1.27  mm  and  stacking  sequence 
[079070790°].  The  material  properties  for  the  composite  plate  are  the  following,  £/=134.4GPa, 
£r=10.3GPa,  v/2=0.33,  GI2=G,3= 5GPa,  G23- 2Gpa  and  y>=1477Kg/m3.  First,  a  single 
delamination  is  introduced  between  the  fourth  and  the  fifth  plies.  The  length  of  delamination  is 
one  third  of  the  length  of  the  plate.  Next,  two  delaminations  are  symmetrically  placed  with 
respect  to  the  plate  edges.  From  Fig.  7,  it  can  be  observed  that  the  proposed  damage  index 
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precisely  indicates  the  delamination  position  and  length,  that  is  elements  with  m= 5-9  and  7=1-4. 
This  damage  index  can  be  used  in  experimental  work  when  the  appropriate  modal  strain 
distribution  is  measured.  Figure  8  presents  the  result  of  the  damage  index  in  the  two- 
delamination  case.  It  can  be  observed  that  the  delaminated  zone  near  the  clamped  end  of  the 
plate  is  precisely  located.  The  elements  with  m= 3,4,5  in  the  global  mesh  have  a  very  high 
damage  index  compared  to  their  neighbors.  For  the  delamination  placed  near  the  free  end, 
elements  with  m- 9,10  are  correctly  indicated  as  belonging  to  the  damaged  zone.  Therefore,  this 
damage  index  is  capable  of  identifying  both  delaminated  zones  quite  accurately. 


in 


Fig.  7  Single  delamination  case 


Fig.  8  Two  delamination  case. 


Next,  the  delaminated  composite  laminate  with  piezoelectric  patches  is  investigated  as 
shown  in  Fig.  9.  The  dimension  of  the  structure  is  of  length  a=240mm,  width  b=135mm  and 
total  thickness  of  composite  zone  h=3mm.  Four  pairs  of  piezoelectric  layers  with  length 
Lp=75mm,  width  lp=45mm  and  thickness  hp=0.5mm  are  assumed  to  be  bonded  at  the  plate 
surface.  The  stacking  sequence  for  this  composite  plate  is  [90°,0o,45o,-45o,0o,90°]s.  The  same 
composite  material  is  used  and  the  piezoelectric  material  properties  are:  £'=63GPa,  i^=0.33, 
G=24.2GPa,  /?=7600Kg/m3  and  J/2=254xl0'12  m/V.  Note  that  no  electrical  force  is  applied  on 
the  piezoelectric  materials,  therefore,  electrical  effects  are  not  considered  in  the  current  study.  In 
the  delaminated  case,  a  single  through  the  width  delamination  of  length  p=45mm  is  assumed  to 
be  placed  between  layers  7  and  8  of  the  composite  plate.  The  numbering  of  layers  starts  with  the 
bottom  layer. 
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Fig.  9  Strain-based  index  (with  delamination  and  actuators) 

As  seen  from  Fig.  9,  the  proposed  strain-based  index  based  is  capable  of  identifying  the 
position  and  extension  of  delamination  quite  accurately.  The  position  corresponds  to  elements 
with  numbers  3,  4  and  5  on  x-axis.  Indeed,  elements  with  numbers  0,  1,  2  and  6-15,  in  the 
longitudinal  direction,  belong  to  either  the  healthy  composite  or  active  undamaged  zones,  while 
elements  with  numbers  3, 4  and  5  belong  to  top  and  bottom  delaminated  sublaminates. 
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Discoveries/Inventions 


The  identification  of  damage  through  the  use  of  dynamic  strain  and  RMS  values  of  the 
strain  is  unique  in  the  field  of  structural  health  monitoring  of  heterogeneous  smart  systems. 
Conventional  health  monitoring  studies  rely  on  expensive  experiments  and  crude  analytical 
models.  The  findings  from  this  research  indicate  that  the  accurate  theory  developed  can  be  a 
useful  tool  in  detecting  the  presence  of  delaminations  in  smart  structures.  Also,  mode  shapes, 
frequencies  and  curvature  mode  shapes  are  commonly  used  in  formulating  damage  indices.  The 
current  research  shows  that  dynamic  strain  is  more  effective  in  capturing  the  location  and  size  of 
delamination  in  composites. 
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